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Abstract 

Within the framework of an isospin and momentum dependent transport model, the emissions 
of isospin particles (nucleons and light clusters) squeezed out in heavy-ion collisions are investi- 
gated as probes of the poorly known symmetry energy at high baryon density. Two different mass 
splittings of neutrons and protons in nuclear medium as m* > m* and m* < m* are used in the 
model and their influence on the isospin emission in heavy-ion collisions is discussed thoroughly. 
The competition between the stiffness and the momentum dependence of the symmetry potential 
on reaction dynamics are compared and systematically analyzed. It is found that the difference 
of the neutron and proton directed flows and the transverse momentum distribution of the neu- 
tron/proton ratio are sensitive to the stiffness of the symmetry energy, which can not be changed 
with the controversial effective mass splitting. The elliptic flows of free nucleons at high transverse 
momentum within mid-rapidity emission are a promising observable as distinguishing the nucleon 
effective mass splitting. 
PACS: 21.65.Ef, 24.10.Lx, 25.75.-q 
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It has been well established the fact that the effective nucleon mass in nuclear matter or 

finite nuclei deviates from its vacuum value [U [2] . Moreover, a splitting of neutron and proton 

effective mass (here, the nonrelativistic mass (Landau mass) is concerned in accordance with the 

definition in Ref. j3j) appears in neutron- rich matter, which increases with the isospin asymmetry. 
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Predictions of the mass splitting based on nuclear many-body theories differ widely. Based on the 
realistic nucleon-nucleon interactions, the Brueckner-Hartree-Fock (BHF) and Dirac-Brueckner- 
Hartree-Fock (DBHF) calculations predict a neutron-proton mass splitting of m* > m* [21 HJ |5] . 
However, the relativistic mean-field (RMF) theory gives a contrary result after the inclusion of 
the isovector mesons p and 5 [HIE]. The Skyrme-Hartree-Fock (SHF) predicts both splittings of 
m* > m* and m* < m* exist with different Skyrme parameters |7j. More realistically, the mass 
splitting of neutron and proton has been constrained from the energy dependence of the Lane 
(symmetry) potential by the nucleon-nucleus scattering experimental data, which moderately 
supports the mass splitting of m* > m* [8j. Further constraints of the effective mass splitting 
from heavy-ion collisions are still necessary and its influence on reaction dynamics would be 
interesting. The Landau mass splitting of neutron and proton in neutron-rich matter results from 
the momentum dependence of the symmetry potential, which directly affects the isospin transport 
in heavy-ion collisions and consequently the extraction of the density dependence of the symmetry 
energy. 

The exact knowledge of the density dependence of the nuclear symmetry energy is still poorly 
known, in particular at high baryon densities, but which has important application in nuclear 
physics itself and also astrophysics. Heavy-ion reactions with neutron-rich beams provide an 
unique opportunity to explore the density- dependent symmetry energy in a broad domain of den- 
sity. A symmetry energy of the form -E'sym(p) ~ 31.6(p/po)'*^ MeV with 7 ^ 0.69 — LOS has 
been roughly extracted for densities between O.lpo and L2po (where po is the nuclear satura- 
tion density) [9J. However, the high-density behavior of nuclear symmetry energy predicted by 
transport models associated with the existing experimental data deviates largely because of the 
inconsistent treatment of the mean-field potentials and the unclear in-medium nuclear interaction 
and properties of resonances [TOl [TTl [121 US], and an opposite conclusion was drawn. Further- 
more, new experimental data related to high-density observables and modifications of transport 
models including the in-medium effects in two-body collisions and in mean-field propagation, are 
very anticipated. In this work, we present systematic investigations of the effective mass splitting 
of neutron and proton and its influence on reaction dynamics. An extraction of the symmetry 
energy from the preequilibrium nucleon emission is performed with an isospin and momentum de- 
pendent transport model (Lanzhou quantum molecular dynamics (LQMD)). The model has been 
successfully applied to treat the nuclear dynamics at near Coulomb barrier energies and also to 
describe the capture of two heavy colliding nuclides to form a superheavy nucleus [131115]. Further 
improvements of the LQMD model have been performed in order to investigate the dynamics of 
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pion and strangeness productions in heavy-ion collisions and to extract the information of isospin 
asymmetric EoS at supra-saturation densities [121 IISI [13 HI [IH] • The momentum dependence of 
the symmetry potential was also included in the model, which results in a splitting of proton and 
neutron effective mass in nuclear medium [20] . 

In the LQMD model, the time evolutions of the baryons (nucleons and resonances (A(1232), 
N*(1440), N*(1535))) and mesons in the system under the self-consistently generated mean-field 
are governed by Hamilton's equations of motion, which read as 

dH . OH 

Pi = -7^, = — . (1) 

Here we only consider the Coulomb interaction for charged hyperons. The Hamiltonian of baryons 
consists of the relativistic energy, the effective interaction potential and the momentum dependent 
interaction. The effective interaction potential is composed of the Coulomb interaction and the 
local potential 

Uint = Ucoul + Uioc- (2) 

The Coulomb interaction potential is written as 

Ucoui = ^ E ^erf{r,,/VAL) (3) 

where the ej is the charged number including protons and charged resonances. The rtj = |rj — rj| 
is the relative distance of two charged particles. 

The local interaction potential is derived directly from the Skyrme energy-density functional 
and expressed as 

Uioc = J Vioc{p{v))dv. (4) 
The local potential energy-density functional reads 



„ „z a n „iso 

vUp) = + V + drP'^VpT + ^(Vp)^ + ^[V(p„ - p,)] 
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where the p„, pp and p = Pn + Pp are the neutron, proton and total densities, respectively, and 
the 6 = {pn — Pp)/{Pn + Pp) is the isospin asymmetry. The coefficients a, (3, 7, g-r, Qsur, gfur 
related to the Skyrme parameters to, ti, ^2, ^3 and Xq, Xi, X2, x^ and the parameter Sly6 is taken 
in the calculation. The bulk parameters a, (3 and 7 are readjusted after inclusion the momentum 
term in order to reproduce the compression modulus of symmetric nuclear matter (here, K=230 
MeV) and the binding energy of isospin symmetric nuclear matter at saturation density, which 



have the values of -296.6 MeV, 197 MeV and 1.143, respectively. The E^ll^ is the local part of the 
symmetry energy, which can be adjusted to mimic predictions of the symmetry energy calculated 
by microscopical or phenomenological many-body theories and has two-type forms as follows: 

Ei"yM = lCsym{p/por% (6) 

and 

Elymip) = asymip/Po) + hsym{p/Po?- (7) 

The parameters Cgym, CLgym and bgym are taken as the values of 52.5 MeV, 43 MeV, -16.75 MeV 
and 23.52 MeV, 32.41 MeV, -20.65 MeV corresponding to the mass splittings of > m* and 
m* < m*, respectively. The values of 7s=0.5, 1., 2. have the soft, linear and hard symmetry energy 
with baryon density, respectively, and the Eq. (7) gives a supersoft symmetry energy, which cover 
the largely uncertain of nuclear symmetry energy, particularly at the supra-saturation densities. 
All cases cross at saturation density with the value of 31.5 MeV. The local part of the symmetry 
energy can be adjusted to reflect the largely uncertain behavior of the symmetry energy at sub- 
and supra-normal densities. 

A Skyrme-type momentum-dependent potential is used in the model as follows: 

Umom^ Y- H ^Cr,r'5r,rA',ri j J J dpdp'drfi{r,p,t) 

x[ln(e(p-p')^ + l)]V.(r,p',t). (8) 
The term is also given from the energy-density functional in nuclear matter, 

Umom = ^ E Cr,r' j j j dpdp'dvU{v, p)[ln(e(p - p')' + 1)] Vr'(r, p')- (9) 

Here Cr,r = Cmomi^ + x), Cry — Cmom{i — x) {j ^ t') and the isospin symbols t(t') represent 
proton or neutron. The sign of x determines different mass splitting of proton and neutron in 
nuclear medium, e.g., positive signs corresponding to the case of m* < m*. The parameters Cmom 
and e was determined by fitting the real part of optical potential as a function of incident energy 
from the proton-nucleus elastic scattering data. In the calculation, we take the values of 1.76 MeV, 
500 c^/GeV^ for the Cmom and e, respectively, which result in the effective mass m* /m—0.75 in 
nuclear medium at saturation density for symmetric nuclear matter. The parameter x is changed 
as the strength of the mass splitting, and the values of -0.65 and 0.65 are respective to the cases 
of m* > m* and m* < m*. 
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Combined Eq. (5) and Eq. (9), we get a density, isospin and momentum-dependent single- 
nucleon potential in nuclear matter as follows: 

UAp, 6, P) = «^ + < + |^?.p^/Vpf + Ei-M5' + ?^^p5^ + Ei-Mp^ 

Pq Pq 6 ■ op OPr 

+ -Cr,r f rfp7.(r,p)[ln(6(p-pr + 1)]^ 
Po -I 

+ -C,y I dp7.'(r,p)[ln(e(p - p')^ + 1)]^. (10) 

Here r ^ t', d6^/dpn = 45pp/p^ and d6^/dpp = —A5pn/p^- The nucleon effective (Landau) mass 
in nuclear matter of isospin asymmetry 5 = {pn — Pp)/{pn + Pp) with p„ and Pp being the neutron 
and proton density, respectively, is calculated through the potential as m* = m^/ (^1 + y^l^l) 
with the free mass rrij. at Fermi momentum p = p^,. Therefore, the nucleon effective mass only 
depends on the momentum-dependent term of the nucleon optical potential. The isovector part 
of the optical potential, i.e., the symmetry or Lane potential [21j, can be evaluated from the 
expression Usym{p,p) = {Un{p,S,p) — Up{p,6,p))/26. Shown in Fig. 1 is a comparison of the 
momentum dependence of symmetry potential with the mass splittings of m* > m* in the left 
window and m* < m* in the right window at the baryon densities of 0.5po, l.Opo and 2.0po, 
respectively. One should note an opposite variation with nucleon momentum appears for the two 
different mass splittings. 

We have included the resonances A(1232), A^*(1440), A^*(1535) in nucleon-nucleon (NN) col- 
lisions and their decaying channels into pions and etas. The in-medium potentials for pion and 
kaon (antikaon) propagations can be chosen in the model p!7l [19] . Furthermore, strange particles 
are directly created by inelastic hadron-hadron collisions [18j. The probability reaching a channel 
in a NN collision is determined by its contribution of the channel cross section to the total cross 
section. The choice of the channel is done randomly by the weight of the probability. 

The preequilibrium nucleons in high-energy heavy-ion collisions are mostly produced during 
a compression stage of two colliding partners within a very short time. Therefore, the high- 
density information of nuclear phase diagram are expected to be extracted from the nucleons 
or light complex particles, which can be constrained from the longitudinal rapidity distributions 
and the azimuthal emissions. Shown in Fig. 2 is a comparison of transverse emission ratios of 
neutron/proton and '^H/^He within the rapidity selection of \y/yproj\ <0.25 in the ^^'^Au+^^^Au 
reaction at the incident energy of 400^ MeV for the near central (b=l fm) and semi-central 
(b=6 fm) collisions, but different mass splitting. One can see that the influence of the symmetry 
energy appears at low transverse momentum and the situation does not be changed with the mass 



splitting. A hard symmetry energy enforces a strongly repulsive force on neutrons in the high- 
density domain, furthermore, squeezes out more neutrons in the preequilibrium stage of dynamical 
evolution. From the negative contribution of the momentum-dependent interaction of m* > m* to 
the potential part of the symmetry energy in Ref. [201, one notices that the momentum- dependent 
potential leads to an attractive force on neutrons, in particular at high densities, which reduces the 
n/p yields at high transverse momentum. Opposite contributions take place and a flat distribution 
of the n/p and ^H/^He appears in the case of m* < m*. The results are nice consistent with the 
stochastic mean-field (SMF) calculations [22J. Experimental measurements of the momentum 
distributions of the n/p ratios would be interesting to get the information of the effective mass 
splitting in nuclear medium and also the symmetry energy at supra-normal densities. 

Collective flow has been verified as a nice approach to reconstruct the reaction plane and to 
study azimuthal correlation of the fireball formed in heavy-ion collisions. Rather systematical 
data has been measured to investigate the equation of state from flow analysis [231 121] • The 
flow information can be expressed as the flrst and second coefficients from the Fourier expan- 
sion of the azimuthal distribution ^{y,Pt) = Nq{1 + 2Vi{y,pt) cos(0) + 2V2{y,pt) cos(20)) [25] . 
where pt = '^Jp'^ + Py and y are the transverse momentum and the longitudinal rapidity along 
the beam direction, respectively. The directed (transverse) flow is deflned as the flrst coefficient 
and expressed as Vi = (px/Pt), which provides the information of the azimuthal anisotropy of the 
transverse emission. The elliptic flow V2 = {{Px/Ptf' ~ (Py/Pty) gives the competition between 
the in-plane (V2 >0) and out-of-plane (V2 <0, squeeze out) emissions. The brackets indicate aver- 
aging over all events in accordance with a speciflc class such as rapidity or transverse momentum 
cut. Shown in Fig. 3 is a comparison of rapidity distributions of directed and elliptic flows of 
free nucleons for near central (b=l fm) and semi-central (b=6 fm) ^^^Au+^^''Au collisions at the 
incident energy of 400^4 MeV with the mass splitting of m* > m* for the hard (left window) 
and supersoft (right window) symmetry energies, respectively. One notices that the difference 
of neutron and proton transverse flows is slight for the near central collisions. The elliptic flow 
disappears for the near central heavy-ion collisions and the influence of symmetry energy on the 
elliptic flow distribution is weak. We calculated the difference of neutron and proton directed 
flows as a function of rapidity distribution in the ^^^Au+^^^Au collisions (b=6 fm) at the energy 
of 400yl MeV with the hard and supersoft symmetry energies, but for the different mass splitting 
as shown in Fig. 4. One can see that an anti-flow spectrum appears in both mass splittings. 
Moreover, the supersoft symmetry energy gives a larger flow difference, in particular around the 
region of projectile (target) rapidity. The phenomena can not be influenced by changing different 



mass splitting. The flows of fast nucleon emissions would be a nice probe of the nuclear symmetry 
energy. 

The effect of mass splitting of neutron and proton on fast nucleon emissions is one of the main 
purposes in this work. Shown in Fig. 5 is a comparison of elliptic flows of free nucleons with 
a complete different mass splitting for the cases of hard and supersoft symmetry energies. One 
notices that the momentum-dependent potential with the m* < m* enhances the out-of-plane 
emissions for neutrons because of its repulsive force on neutrons. Thus, the flow is reduced and 
even coincides with protons. Furthermore, a larger mass splitting (stronger isovector momentum- 
dependent interaction) increases the squeeze out neutrons and consequently results in an lower 
elliptic flows for neutrons [22j. The stiffness of the symmetry energy slightly affect the flow 
spectrum. Again, we compared the transverse momentum distribution of elliptic flows of proton 
and neutron as shown in Fig. 6 with the mid-rapidity selection of \y/yproj\ <0.25. A similar effect 
of the effective mass splitting with the rapidity distribution is observed. Basically, the symmetry 
energy does not change the spectrum structure. The elliptic flows of hydrogen-like particles have 
been used to extract the high-density symmetry energy [13j. More sensitive observable can be 
seen from the flow difference between neutron and proton as shown in Fig. 7. The m* > m* 
case is always larger, in particular in the domain of high transverse momentum. The findings are 
nice consistent with the SMF calculations [22] although different treatment of the momentum- 
dependent potential. 

In summary, within the transport model (LQMD) we have investigated the effect of the 
momentum-dependent interaction on fast nucleon emission. Two different mass splittings, i.e., 
m* > m* and m* < m*, are used in the calculations. The momentum-dependent potential plays 
a significant role on the fast nucleon emissions. Specifically, the high-momentum neutron-proton 
ratio is changed, which appears more sensitive to the mass splitting than to the stiffness of the 
symmetry energy. The neutron/proton ratio and the directed flow difference of neutron and pro- 
ton via transverse momentum distribution in the domain of mid-rapidity are nice observables for 
extracting the high-density symmetry energy. 
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Figure 1: Momentum dependence of symmetry potential at baryon densities of 0.5po, 1-Opo and 
2.0po with the mass sphttings of m* > m* (left panel) and m* < m* (right panel), respectively. 
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Figure 2: Transverse momentum distributions of neutron/proton and ^H/'^He ratios within the 
rapidity bin \y/yproj\ <0.25 in the ^^''Au+^^^Au reaction at the incident energy of 400 MeV/nucleon 
for the near central (b=l fm) and semi-central (b=6 fm) collisions with different mass splitting. 
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Figure 3: Rapidity distributions of directed and elliptic flows of free nucleons in the Au+^^^Au 
reaction at the incident energy of 400 MeV/nucleon for the near central (b=l fm) and semi-central 
(b=6 fm) collisions with the mass splitting of m* > m* for the hard (left panel) and supersoft 
(right panel) symmetry energies, respectively. 



19 



0.10 



0.05 



0.00 



-0.05 



-0.10. 




proj 



proj 



Figure 4: Comparison of the difference between neutron and proton directed flows in semi-central 
^^^Au+^^^Au collisions at the energy of 400 MeV/nucleon with the hard and supersoft symmetry 
energies for the different mass splittings of m* > m* (left panel) and m* < m* (right panel), 
respectively. 
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Figure 5: Rapidity distributions of elliptic flows of free nucleons with the different mass splittings 
at the cases of hard (left panel) and supersoft (right panel) symmetry energies. 
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Figure 6: The same as in Fig. 5, but for transverse momentum distributions with the rapidity bin 

\y/yproj\ <o.25. 




Figure 7: Comparison of the difference between neutron and proton elliptic flows with the rapidity 
bin \y/yproj\ <0.25 for the different mass splitting in the semi-central ^^'''Au+^^'^'Au coUisions. 
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